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We derive new limits on light stops from diboson resonance searches in the 77 , 
Zj, ZZ, WW and hh channels from the hrst run of the LHC. If the two-body decays 
of the light stop are mildly suppressed or kinematically forbidden, stoponium bound 
states will form in pp collisions and subsequently decay via the pair annihilation of 
the constituent stops to diboson hnal states, yielding striking resonance signatures. 
Remarkably, we hnd that stoponium searches are highly complementary to direct 
collider searches and indirect probes of light stops such as Higgs coupling measure¬ 
ments. Using an empirical quarkonia potential model and including the hrst two 
5-wave stoponium states, we hnd that in the decoupling limit < 130 GeV is 
excluded for any value of the stop mixing angle and heavy stop mass by the combi¬ 
nation of the latest resonance searches and the indirect constraints. The 77 searches 
are the most complementary to the indirect constraints, probing the stop “blind 
spot” parameter region in which the trilinear coupling is small. Interestingly, 

we also hnd that the Z'j searches give a stronger constraint, < 170 GeV, if the 
stop is primarily left-handed. For a scenario with a bino LSP and stop NLSP, several 
gaps in the direct collider searches for stops can unambiguously be hhed with the 
next run of the LHC. For a stop LSP decaying through an R-parity violating UDD 
coupling, the stoponium searches can hll the gap 100 GeV < mi^ < 200 GeV in the 
direct searches for couplings A" < 10“^. 
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I. INTRODUCTION 

Light stops are a quintessential feature of a natural supersymmetric theory HHn, being 
responsible for the cancellation of the dominant Higgs mass quadratic divergence coming 
from the top quark. However, after Run 1 of the LHC the possible existence of light stops 
has been strongly constrained by a suite of dedicated searches by ATLAS and CMS [SHE]. 
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For instance, in a simplified scenario containing a nentralino as the lightest snpersymmetric 
particle (LSP) and the stop as the next-to-lightest snpersymmetric particle (NLSP), with 
a few gaps and caveats, stop masses as high as ~ 700 GeV and LSP masses as high as 
~ 250 GeV have been excluded mn- Still, as is always the case with direct searches, 
the limits depend strongly on the spectrum and decays of the stops, and there remain open 
windows in which a light stop can hide from LHG searches. Examples include compressed [H] 
or stealth [10] stops and R-parity violating stop decays mm- It is therefore critical to 
continue to develop new strategies to directly search for light stops 0 . 

Another opportunity to probe light stops is presented when the stop is long lived, which 
naturally occurs in a number of motivated scenarios. Indeed, provided the stop does not 
have an unsuppressed 2 -body decay, stop pairs produced through gluon fusion can form a 
stoponium bound state due to the Goulombic attraction mediated by the strong force 
\n\ . Once the bound state is formed, it can decay via annihilation to (y' 5 ', 77 , Z-j, WW, ZZ 
and hh, leaving a distinctive resonance signature. As we will emphasize below, the branching 
ratios to these hnal states depend essentially on one parameter, the trilinear Higgs-stop-stop 
{h^titl) coupling. Therefore, provided there is enough time for the stoponium to form and 
annihilate decay, the signature is independent of the details of the light SUSY spectrum and 
thus offers a highly complementary probe to direct stop searches. 

In this paper we will use the null results from 77 , Z-j, WW, ZZ and hh resonance searches 
to constrain the production of stoponium during the hrst run of the LHG and in turn 
derive new limits on light stops. Assuming the stop is sufficiently long-lived, constraints 
on stoponium production in the 77 , 7 Z, and ZZ channels alone limit stops lighter than 
mi^ < 125 GeV for an^ value of the stop mixing angle and heavy stop mass mi^. Furthermore, 
for mostly left-handed light stops, the Z 7 searches constrain light stops up to about 170 
GeV. We will also present projections for 14 TeV LHG with low (30 fb~^) and high (3 ab“^) 
luminosities. 

One new aspect of our study which has not been emphasized in the literature is the 
interplay between stoponium searches and other indirect constraints on light stops. In 
particular, since the stoponium branching ratios depend primarily on the h^tit\ coupling, 
our new constraints are highly complementary to those suggested by Higgs coupling mea¬ 
surements. Light stops with a large coupling to the Higgs boson signihcantly alter Higgs 
production in the gluon-gluon fusion channel and are already strongly constrained except 
for a “blind spot” in the parameter region where the h^tit\ coupling is small [2^ 123] . On 
the other hand, the 77 resonance signature of stoponium is enhanced in the region where 
the h^tit\ coupling is small. Gombining both constraints, we can rule out the existence of 
light stops up to about 130 GeV independently of the other parameters in the stop sector. 
Besides the interplay with Higgs couplings, we will also explore the impact of other indirect 


4 


constraints on light stops, inclnding vacunm stability and electroweak precision data. 

We also discnss the conseqnences of onr resnlts for two motivated scenarios in which sto- 
poninm annihilation decays are relevant. The hrst is the canonical R-parity conserving case 
with a bino LSP and a stop NLSP. For mi^ '^rnt + m^o the stoponinm annihilation decays 
are visible, providing a clean and nnambignous probe of light stops which is independent 
of possible kinematic degeneracies or model dependencies of stop branching ratios. Second, 
we discnss the case of a stop LSP decaying via the R-parity violating UDD operator to 
two jets, which is challenging to probe directly dne to the large QCD backgronnds. Again, 
stoponinm can provide a clear-cnt test of light stops in this scenario. 

It is worth emphasizing that there is a potentially large sonrce of theoretical nncertainty 
present in the stoponinm prodnction cross section coming from the assnmed potential model 
and contribntion of the excited states to the signal. The limits we present here are based on 
the empirical qnarkonia potential model of Hagiwara et al. [21]. Fnrthermore, as discnssed in 
Ref. [20|, it is likely that higher F-wave stoponinm states contribnte to the resonance signals, 
and we therefore inclnde the hrst two states in onr estimate. As we will argne below, we 
believe onr assnmptions in this regard are conservative, bnt we will also discnss in detail 
how different choices affect the stoponinm prodnction cross section and the resnlting limits. 

We begin in Section [IT] with a brief review on the prodnction and decay of stoponinm at 
hadron colliders. Onr new limits and fntnre projections from LHC resonance searches on 


several benchmark scenarios for each diboson channel are presented in Sec. Ill The cnrrent 
limits from stoponinm and other indirect constraints on the fnll stop sector parameter space 


are combined in Sec. IV In Sec. Mwe discnss the implications of onr results on two motivated 


scenarios containing lights stops. Section VI presents our conclusions and outlook. Finally 
we include several Appendices which collect and summarize our conventions for the 
stop sector, the indirect constraints on light stops, the stop decay widths, our procedure for 
limit extrapolation, and the LHC diboson resonance searches used in our analysis. 


II. STOPONIUM PRELIMINARIES 

In this work we will investigate the annihilation decay signatures of the S'-wave {J^^' = 
stoponinm bound state, r]i, at the LHC. These signatures result from the following 
processes: 

pp^r]-t + X, Pi ^ 77, Z7, VFIV, ZZ,/i/i... ( 1 ) 

Searches for stoponinm resonances provide a highly complementary probe of light stops, 
testing different assumptions about SUSY spectra and parameters than both direct searches 
for stops and indirect probes such as Higgs coupling measurements, precision electroweak 
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data, and vacuum stability constraints. In particular, direct searches for stops (and other 
new particles) are inherently model dependent, requiring assumptions about the superpart¬ 
ner spectrum and stop decay chains. In contrast, the limits and projections we will derive 
for stoponium will depend on two different assumptions: 1) the stop does not have an un¬ 
suppressed 2-body decay, and 2) the value of the trilinear Higgs-stop-stop coupling. As we 
will discuss shortly, the hrst assumption of a narrow stop width is actually satished in a 
number of interesting scenarios. The second assumption on the value of the h^titi coupling, 
which is determined by the stop sector masses and mixing angle, will govern the stoponium 
branching ratios for the various hnal states in Eq. ([^. 

In this section we will collect the relevant ingredients entering into our study of stoponium. 
Our analysis relies on a number of previous works [HHIZ!, notably those of Martin and 
collaborators [IHH2I]. Next we describe in more detail the basic requirements for observable 
stoponium annihilation decays. Following this we will review the production of stoponium 
at the LHC and the stoponium decay branching ratios. 

Before beginning, let us note here that we will work exclusively in this paper with the 
physical stop parameters Ot}, in which ii (^ 2 ) denotes the lightest (heaviest) phys¬ 

ical stop state, and a mixing angle 9 = 0 (vr/2) corresponds to a purely left-handed (right- 
handed) stop. We will also restrict to the decoupling limit. We refer the reader to Ap¬ 
pendix 1^ for our stop sector conventions. 


A. Conditions for stoponinm formation and annihilation decays 

What conditions are required for the distinctive resonance signatures of stoponium to 
be observable? ^ Provided they are kinematically accessible, stop pairs will be copiously 
produced in pp collisions. However, whether or not the bound state forms and subsequently 
decays via the annihilation processes in Eq. ([^ depends on the natural width of the stop. The 
bound state formation time is on the order of the inverse binding energy which, assuming a 
Coulombic potential, is given by Eb ~ as{aQ^Y ~ 0(1 GeV), where oq ~ [as{aQ^)mi]~^ 
is the “Bohr radius” of the system. Therefore, stoponium will form provided the natural 
width of the stop is smaller than the binding energy. Subsequently, the decay of the bound 
state may proceed either through the prompt decay of the constituent stop or through the 
annihilation decays. For instance, if the coupling is small the dominant annihilation 


^ Useful related discussions for toponium bound states can be found in Refs. nans]. 
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decay is to a pair of gluons, with partial width (see also Sec. II C| below) 


^{Vt^ 99) - o^s- - 


m 


Vi 


( 2 ) 


where i?(0) = \/^'ijj{0) is the radial wavefunction at the origin, for which we employ the 
parameterization in Ref. [IH1I21]- For mi^ ~ 0(100) GeV, the factor |R(0)p/m^, ~ 0.2 GeV, 
leading to a partial width T{t]i —)• gg) ~ 0(1) MeV in the light stop mass range of interest. 
If the annihilation decay width dominates over the natural stop width we can potentially 
observe the stoponium through its resonance signature. 

In summary we require, 

Fij ^ F^-,, Eh. (3) 

In practice, since the annihilation decay width is typically smaller than the binding energy, 
it sets the upper bound in Eq. (|^. In fact, these conditions are satisfied in a number of 
interesting of scenarios and over a wide range of model parameters, particularly when the 
stops are light. We will examine two such scenarios in Section fV] 


B. Production 


At leading order (LO) the production of stoponium in pp collisions at the LHG proceeds 
via gluon-gluon fusion, with a cross section given by [IH] 


vr 


(^Lo{PP Vi) = ^^^iVi 99)'Pgg{'r), 


8m3 


(4) 


where Vgg is the gluon parton luminosity as a function of r = rrig-/s, with s the squared 
center-of-mass energy, dehned in Eq. (D.2) in Appendix]^ The partial decay width for 
F(^t —t gg) is given in Eq. (|^. Next-to-leading order (NLO) perturbative QGD corrections 
to stoponium production have been computed in Ref. m and with the threshold resum¬ 
mation in Ref. |T7]. To account for these corrections we will use the results of Ref. [20]. In 
particular, we extract the cross section prediction from their Fig. 9 and furthermore estimate 
a conservative ~ 10% scale uncertainty on this prediction from their Fig. 7 |20|. We can also 
reasonably estimate a ~ 5 — 10% PDF uncertainty in the mass range of interest, similarly 
to that for Higgs production via gluon-gluon fusion [26] . 

There are likely larger sources of theoretical uncertainties than the QGD scale uncertainty, 
which are, however, difficult to quantify. Notice in particular that the cross section in Eq. (|^ 
depends on the value of the stoponium wavefunction at the origin squared, |R(0)p, through 
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Eq. ([^, and our imprecise understanding of this matrix element potentially provides a 
significant source of uncertainty. Ref. [2D] employs the = 300 MeV parameterization of 
the wavefunction at the origin from the study of Hagiwara et ah [2l|. This parameterization 
is based on an empirical quarkonia potential model, which deviates signihcantly from the 
pure Coulombic form at large due to the effect of higher order QCD corrections. For 
instance, the prediction for |i?( 0 )p for the ground state is smaller by roughly a factor of 2 
in comparison to that obtained from a Coulomb potential, suggesting that our estimates of 
the production rate are conservative. Ref. |2l] also present = 200 MeV and A^ = 400 
MeV parameterizations, which lead to a ~ 10% reduction or enhancement in the stoponium 
production rate, respectively, giving some sense of the uncertainties involved. We note that 
the study of Ref. m is 25 years old, and given the signihcant experimental and theoretical 
progress in the studies of quarkonia and QCD in the interim, it would be extremely useful 
to perform a modern analysis in the spirit of [ 21 ], with the particular aim of assessing the 
theoretical uncertainties. 

Another source of uncertainty in the prediction is the contribution of higher R-wave 
stoponium states to the annihilation decay signals. Ref. [2D] includes the hrst two states in 
their calculation, but reasonably argues that higher states will likely add to the signal, and in 
this way our estimates of the production rate are likely conservative. In Figure we display 
the stoponium production cross section at the LHC under several different assumptions 
regarding the potential and contribution of excited states. Needless to say, it would be 
extremely interesting to study the stoponium system in more detail to both quantify the 
uncertainty due to the modelling of the bound state potential as well as understand the 
phenomenological consequences of the higher radial states. 


C. Decay branching ratios 

For the stoponium decay branching ratios, we use the LO results presented in Ref. [18] 
for all decay channels except for the 77 ^ —)■ 77 channel. For the diphoton channel we obtain 
the partial width T{r]f —)■ 77 ) by multiplying the LO partial width ^LoiVi 99) [IB] by 
the NLO ratio rNLo(ht 77)/rNLo(ht gg) computed in Ref. [19]. This approximation 
reproduces the correct NLO BR( 77 ) when BR( 5 f 5 f) dominates, and it is this case in which 
the diphoton channel is most important^. When other decay modes such as hh or WW 
dominate instead, the diphoton branching ratio is anyway too small to induce an observable 
signal. We use the two-loop running as{g) and the one-loop running a{g) renormalized 

^ Although the Z 7 channel is also important in this parameter region, the NLO corrections for this channel 
are not available. 




FIG. 1: Stoponium production cross sections at the LHC. Our baseline cross section for ^/s = 8 
TeV (solid blue) and = 14 TeV (solid green) from Ref. [20], which assume the = 300 
MeV parameterization of the stoponium wavefunction at the origin from Ref. |23| and include 
the contribution of the first two 5-wave states. These cross sections are used to derive our limits 
and projections. For ^/s = 8 TeV we also display for comparison the cross section derived with 
the quarkonium potential | 23 | including instead only the ground state (blue dashed) and the first 
10 5-wave states (blue dotted), as well as the cross section assuming a Coulomb potential (red) 
(including only the ground state) with as evaluated at /r = l/{ris) = 2 a 5 (/u) . 


from as{mz) = 0.118 and a{mz) = 1/128.0. 

We will work in the deconpling limit. The minimal low energy spectrnm will contain the 
stops and possibly the left-handed sbottom. It is also possible that the spectrnm contains a 
bino-like nentralino LSP, since the condition ^ is readily satished when < rrit+m^o (see 
Sec. V for a detailed discnssion of this scenario), althongh the stoponinm annihilation decays 
to nentralinos in this case are always negligible. We will not consider Higgsinos or Winos that 
are lighter than ti, as in this case we expect the stop to have an nnsnppressed two body decay 
to a lighter chargino and a bottom qnark, generically violating the condition All other 
heavy Higgs scalars and snperpartners are taken heavy and deconpled and will not inflnence 
the stoponinm decays. The mass of the left-handed sbottom is hxed by assnming vanishing 
sbottom mixing (see Eq. (B.8)). We emphasize that left-handed sbottoms can have non- 
negligible effects on the stoponinm branching ratios when the lightest stop is primarily left- 
handed. As typically happens, nnitarity and gange invariance indnce a (partial) cancellation 


^ In any case, such scenarios with the Higgsino or Wino lighter than the stop are already strongly constrained 
by searches of bb plus missing energy m- 
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FIG. 2: The coupling )^hoi^i*/v (upper) and stoponium branching ratios (lower) for = 120 GeV 
(left) and 160 GeV (right) as functions of the stop mixing angle. A strong correlation between 
and the stoponium branching ratios is observed. The shaded regions are constrained indi¬ 
rectly by Higgs signal strength data, precision electroweak measurements, and vacuum stability as 
discussed in Appendix Here we have fixed = 1 TeV. Note that in the left panel we have 
terminated the plot at small 6t where the sbottom becomes tachyonic. 


in the WW partial width between t-channel stop and sbottom exchanges if the stop is mostly 
left-handed na. Finally, the stoponium decays to bb and tt are also subdominant and will 
not be considered further. 

The main decay modes, rji — )■ gg, 77 , Z 7 , WW, ZZ and hh, are described well by the 
three physical stop parameters, 6 t, and = tan/3. These parameters are defined 

in Appendix The dependence on is weak in most cases, so we £x tp = 10 throughout 
in this paper. The trilinear coupling of the 125 GeV Higgs to the light stops, plays a 

crucial role in determining the stoponium decay pattern. This coupling is written in terms 
of the physical stop parameters as 




V2v 


mlc2p 

v2 ^2 




+ 


mj 

, 2„2 ti 



(5) 
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Due to the last term in Eq. ([^, this coupling can become much larger than the electroweak 
scale, V = 174 GeV, when the stop mixing is large, inducing sizable partial decay widths for 
stoponium to WW, ZZ, hh via the s-channel exchange of h^. This is illustrated in Fig. 
where we display the mixing angle dependence of the coupling Xhoi^i* and the stoponium 
branching ratios. A strong correlation is observed between the value of and the 

pattern of branching ratios. For large values of mixing, the coupling grows and the 

stoponium dominantly decays to WW, ZZ, and, if kinematically allowed, hh. Instead, when 
Xho^j* is small, the gg decay dominates and, importantly, the 77 and Z 7 branching ratios 
are maximized (see Sec. Ill A and Sec. IIIB for further discussion). 

Large values of are constrained, particularly for light stops, by indirect tests includ¬ 

ing Higgs signal strength measurements, electroweak precision data, and vacuum stability. 
These indirect constraints are reviewed in Appendix and shown as the shaded regions in 
Fig. In particular, we observe that only a small region around = 0 is consistent 

with these indirect constraints for very light stops. It is therefore very interesting to consider 
the rji —?• 77 (and Z 7 ) channel as it can provide an independent probe of this unconstrained 
region of parameter space. 


III. NEW LIMITS ON LIGHT STOPS AND FUTURE PROSPECTS 


We now derive the current limits on light stops from stoponium annihilation decays and 
estimate the future reach of each diboson channel in Eq. ([^: 77 , Z 7 , hhlU, ZZ and hh. 
In this section, we simply assume that the stoponium forms and decays via annihilation, 
but we will discuss particular models and relevant parameter space in Sec. |V| For each 
diboson channel, we dehne and consider the ideal or benchmark branching ratio that best 
represents the optimal or characteristic reach, respectively, over a wide range of the stop 
masses. Based on these benchmarks, we present the limits as a function of the stoponium 


mass. Next, in Sec. IV we exhibit the exclusion limits in the general stop parameter space, 
he., as a function of the stop masses and mixing angles, accounting for the realistic branching 
ratios and comparing with the indirect constraints described in Appendix 

Our limits and projections are obtained by recasting the results of the latest heavy res¬ 
onance searches at LHC 7-1-8 TeV, which are collected in Appendix]^ We will present the 
strongest limit in each diboson channel. For the future sensitivities, we extrapolate the 
current expected limits to LHC 14 TeV with 30 fb“^ and 3 ab“^. Our extrapolation method 
and assumptions on statistical and systematic uncertainties are described in Appendix 
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FIG. 3: Current bounds (upper) and future sensitivity (lower) on the ideal 77 scenario defined in 


Sec. Ill A ATLAS (right) and CMS (left) 77 resonance search results are used; the CMS result is 


extended to a higher mass. Four theoretical predictions are shown as in Fig. considering the 
quarkonium-inspired potential of Ref. [23] and including first two 5'-wave states (thick solid), IS 
only (bottom dashed), first 10 5-wave states (mid-dashed) as well as considering a IS with Coulomb 
potential (top dot-dashed). A ±15% scale±PDF uncertainty is also shown as a band around the 
thick line. Vertical shaded regions are excluded if 1S±2S states are considered. A constant fiducial 
cut efficiency 0.66 is unfolded for the ATLAS results. 


A. 


77 


The diphoton channel has long been recognized as a promising place to search for stopo- 
nium [15]. At LO, the —)■ gg and gi —)■ 77 partial decay widths do not depend on stop 
sector parameters other than the lightest stop mass, They are induced solely by the 
strong and electromagnetic interactions, such that the ratio. 


r(77) 


9«r 


( 6 ) 


is fixed by gauge couplings and independent of . As discussed in Sec. |II C 


m our numerics 


we 


include the NLO corrections to the ratio in Eq. (g from Ref. [IHlEn]. These corrections, 
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along with the running gauge couplings, induce a mild dependence on the lightest stop mass, 
with the ratio varying from 0.0025 — 0.0042 as the lightest stop mass varies from 
100 - 400 GeV. The NLO corrections actually reduce the ratio 0 by 35%-25% from the 
LO result pT] . 

Since the gg and 77 decay modes are always present and the size of the corresponding 
partial decay widths are rather model-independent, the diphoton branching ratio is maxi¬ 
mized simply when the other possible decay channels (WW, ZZ, hh) are minimized. We 
will therefore dehne the ideal 77 scenario by the idealized limit in which the branching ratio 
BR( 7 ^ —)■ 77 ) is given by the ratio in Eq. (|^ (including NLO corrections as discussed above). 
We note that the idealized branching ratios are 10^ - 10® times bigger than that of the SM 
Higgs boson with the same mass in the range of = 200 — 600 GeV. 

In the full stop parameter space, the ideal diphoton scenario is approximately realized 
in the region where the Higgs-stop-stop coupling is small, corresponding to mixing 

angles 



In this region of parameter space the —)■ WW, ZZ, hh diagrams involving s-channel Higgs 

exchange are of a similar size or smaller than those mediated by the weak gauge couplings, 
thus explaining the dominance of the gi —>■ gg mode. This can also be seen from the 
dependence of A^joy^y* and the stoponium branching ratios on the stop mixing angle in Fig. 
where we observe that the gg branching ratio is nearly 100% when the is small. For 

the stop mass = 100 — 300 GeV = 1 TeV), the ideal 77 scenario is approximately 
realized for 6 t = 0.10 — 0.22 or 6 t = 1.32 — 1.48. The actual maximum diphoton branching 
ratio in these parameter regions is smaller than the ideal one of Eq. (§ only by about 
0.7 —1.8% fraction. In the next section, we will translate our limits to the full stop parameter 
space, accounting for the realistic branching ratios. 

We note that the ideal 77 scenario lies within the blind spot region consistent with 
current indirect constraints, even for very light stop masses of order 100 GeV. In particular, 
the constraints from the Higgs signal strength data are relaxed in this region because the 
coupling A^joy^y* is small and the modihcation of the gluon-gluon fusion cross section from 
stop loops is minimized. This makes the gi —)■ 77 channel particularly interesting to study 
as it can provide a complementary probe of this parameter region. 

In Fig. 1^ we show the current limits from ATLAS and GMS diphoton resonance searches 
and the future sensitivities for the ideal diphoton scenario. Stoponium masses in 

the the vertical blue bands are excluded by the observed limits. We note that the ATLAS and 
GMS observed limits display complementary patterns of statistical fluctuations, excluding 
somewhat different stop mass ranges. The union of the two searches excludes stop masses up 
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to about ~ 136 GeV. We emphasize that these exclusions are based on the assumption 
of including the hrst two S'-wave stoponium states in the production rate, as described in 
Sec. |II B This is a conservative assumption as it is likely that higher modes also contribute 
to the diphoton signal. Following Ref. P!. for comparison we have also displayed in Fig. 
the cross section including only hrst S'-wave mode or hrst 10 S'-wave modes, allowing the 
reader to see the ehect of this uncertainty. The reader can also see the diherence between 
the two potential models: the thick blue line indicates the empirical quarkonia potential in 
Ref. [21] and the dot-dashed red line represents the Coulomb potential. 

We also show in Fig. [^the 14 TeV projections, based on the extrapolation of the current 
expected limits; see Appendx for a description of the methodology behind our extrapo¬ 
lation. At the early stage of running with 30/fb, stops as heavy as about 180 GeV can be 
excluded in the ideal diphoton scenario, while a high luminosity run with 3/ab can exclude 
the stop up to about 370 GeV. 


B. Zy 


Another important channel is Zy since the branching ratio is also maximized when the 
coupling XfiOi^q is small and the indirect constraints are absent. As for the diphoton channel, 
we can write the ratio of Zy to gg partial widths at LO: 


r(^ 7 ) _ ^ 1 ( 2 _ 4 2 y / ^ ^ 

ngg) ~ r 3 


( 8 ) 


As in the diphoton case, the BR(Zy) is maximized when the BR( 5 f 5 f) dominates and the 
other decay modes {ZZ,WW,hh) are suppressed, corresponding to small Xhohq- However, 
unlike the diphoton case the mixing angle dependence in Eq. ([^ implies that the BR(Zy) 
is largest when the lightest stop ti is mostly left-handed. This asymmetry is also clearly 
visible in Fig. We therefore dehne the ideal Zy scenario by evaluating BR(Zy) with stop 
parameters corresponding to Xf^oqq = 0 and a mostly left-handed ii. In practice, we choose 
different mixing angles satisfying this condition for different stop masses. The mixing angle 
and the BR(Zy) for the ideal Zy scenario gradually increase from 6 t = 0.169 to 0.184 and 
from 0.84% to 1.37%, respectively, as the stop mass varies in the range mq = 100 — 400 
GeV. We note that these branching ratios are 100 - 3000 times larger than those of the SM 
Higgs boson with the same mass in the range of = 200 — 600 GeV. 

We show the current limits on the ideal Zy scenario from the Zy —)■ f+f^y resonance 
searches in the upper panels of Fig. Notably, the GMS and ATLAS results differ somewhat 
in the excluded stop mass range; the GMS result is stronger and excludes stops up to about 
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FIG. 4: Current bounds (upper) and future sensitivities (lower) on the ideal Z 7 scenario defined in 
Sec. IIIB from ATLAS (right) and CMS (left) results. Theoretical predictions and excluded regions 


are as in 


Fig-i A constant fiducial cut efficiency 0.6 is unfolded for ATLAS results. 


170 GeV in the ideal Z 7 scenario, while instead the ATLAS result is sensitive only up to 
120 GeV stops. In terms of the excluded cross-section the GMS result is about 2-3 times 
stronger for the stop mass range shown. Since the Z 7 result can potentially give the strongest 
constraint in the currently unconstrained parameter space with a mostly left-handed stop 
and small \hHiii i we encourage further dedicated analysis between the two experiments. We 
also show the 14 TeV prospects in the lower panels of Fig. The GMS-extrapolated result 
can potentially probe the stop up to 225 GeV with 30/fb and the ATLAS-extrapolated result 
can be sensitive to 300 GeV with 3/ab. As we will see in Fig. in the next section, the 
diphoton channel is most complementary to the indirect constraints, as it is sensitive to the 
entire small region including primarily right-handed stops. On the other hand, the 

Z 7 channel provides stronger limits than the diphoton channel for mostly left-handed stops. 

It is important to mention that our Z 7 results are based on the LO calculation while the 
77 result in the previous subsection is from the NLO calculation in Ref. [ISlEo]. For the 
latter case, the NLO effects reduce the ideal 77 branching ratio by 35%-25% for the stop 
mass in the range 100-400 GeV. If the same NLO reduction is applied to the ratio of the 
Z 7 to gg partial widths (hence, to the ideal BR(Z 7 )), the expected exclusion on the stop 
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FIG. 5: Current(left) and future(right) bounds on the benchmark ZZ scenario with BR(ZZ)=0.06. 
CMS ZZ ^ 4i search result is used. 
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FIG. 6: Current(left) and future(right) bounds on the benchmark WW scenario with 
BR(VFVF)=0.12. CMS WW —)• 2i2iy search result is used. 


mass is approximately reduced by about 20 GeV. In addition to the NLO perturbative QCD 
effects, uncertainties from the number of S'-wave stoponium states contributing to the signal 
and from potential models are again significant as shown by multiple theoretical predictions 

in Fig. 1^ 


ZZ and WW 


The branching ratios BR(ZZ) and BRiWW) are intimately related with the coupling 
as discussed in Sec. IIC The indirect probes of light stops discussed in Sec. [TVj place 
strong constraints on large values of this coupling and in turn determine the maximum 
phenomenologically allowed branching ratios BR(ZZ) and BR(lFfF), as is depicted in Fig.|^ 
Unlike the diphoton channel discussed above, there are several reasonable ways to de¬ 
fine the benchmark ZZ^ WW scenarios. For instance, one could simply take the largest 
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theoretically allowed branching ratio, or alternatively the largest branching ratio consistent 
with indirect constraints; in both cases one is led to vary the mixing angle as a function of 
the stop masses. One could also choose a hxed mixing angle for a range of the stop mass. 
For simplicity, we will instead choose a prescription with a hxed branching ratio that is 
broadly consistent with the current indirect constraints. Referring to Fig. |^of = 160 
GeV case, we dehne the benchmark ZZ, WW scenarios by choosing BR{ZZ)=6% and 
BR(hFhF)=12%. In fact, for the range of 6 t allowed for the 160 GeV stop we obtain maxi¬ 
mum BR(ZZ) ~ 4 - 8% and BR{WW) ~ 11 - 16% for the stop mass = 100 — 300 GeV 
(with = 1 TeV). Thus, the maximum branching ratios in this range of dt do not vary 
much with the stop mass. The current indirect constraints are somewhat weaker for heavier 
stops and a wider range of dt is viable, but we envisage that indirect constraints will also 
improve with more data, probing a broader range of mixing angles and heavier stop mass. 
We refer to the = 160 GeV case as a representative benchmark scenario in view of the 
interplay between stoponium and indirect constraints. Again, we provide limits on the full 
stop parameter space accounting for true branching ratios in the next section. 

Using a hxed branching ratio as a benchmark scenario for a wide range of 0(100) GeV 
stop masses is also consistent with the softened Goldstone enhancement of the stop decay 
into longitudinal Z and W bosons, which could have sensitively increased the branching 
ratio with the stop mass. The decay into gauge bosons are maximized when the amplitude 
is dominated by the s-channel exchange of Higgs boson, corresponding to values of the 
coupling strong enough to dominate the decay but weak enough to be consistent with indirect 
constraints. But this contribution softens the Goldstone enhancement factor rn^z) of 
the amplitude by a suppression l/(m^_ — m^) ~ l/(4m|) from the propagator, leaving the 
dominant amplitude relatively insensitive to the stop mass for the benchmark parameter 
space. 

In Fig. 1^ and Fig. we show the current limits and prospects of the benchmark ZZ 
and WW scenarios from the ZZ ^ M and the WW —?• 2£2z/ channels, respectively. The 
constraint from the ZZ channel is stronger, just as it is for the SM Higgs. The WW —>■ iv2j 
channel gives similar bounds as the 2£2z/ channel but is not shown for simplicity. Gurrently, 
the ZZ channel constrains stops up to about 110 GeV and a small range of 121-126 GeV 
for the benchmark scenario. In the future, the ZZ channel is expected to probe 140 (250) 
GeV stops at i/s = 14 TeV with 30/fb (3/ab) for the ZZ benchmark scenario. While there 
are no limits yet, the LHG at 14 TeV with 3/ab is expected to probe 180 GeV stops in 
the benchmark WW scenario. We note that in the full stop parameter space, with realistic 
branching ratios, there are in fact new limits from the WW channel for very light stops, 
where the hh mode is absent. 


17 



250 300 350 400 450 500 


^stoponium“2/71stop [GcV] 



250 300 350 400 450 500 


^stoponium“2/71stop [GcV] 


FIG. 7: Current(left) and future(right) bounds on the benchmark hh scenario with BR(/i/i)=0.25. 
ATLAS hh —)• 6677 search result is used. 


D. hh 

Similarly to the ZZ and WW branching ratios, the maximnm branching ratio Br(hh) is 
limited by indirect constraints on the conpling Xhoij*. We define the benchmark hh scenario 
by a fixed BR(hh)=0.25. The actnal maximal BR(hh) for the range of 9t cnrrently allowed 
for the 160 GeV stop gradually increases from 16% to 41% for the 130 GeV to 300 GeV 
stop. 

In Fig. we display the current sensitivity and 14 TeV projections on the benchmark 
hh scenario from the hh —)■ bb'j'j resonance search. No new limits can be obtained with this 
channel from the Run 1 data. However, in the long term this channel is expected to probe 
stops as heavy as 215 GeV in the benchmark hh scenario at ^/s = 14 TeV with 3/ab. 


IV. SUMMARY: INTERPLAY OF STOPONIUM AND INDIRECT 

CONSTRAINTS 


Our main results are presented in Fig. which displays the constraints from diboson 
resonance searches on stoponium in the physical stop parameter space. We show six slices 
of the parameter space, fixing the lightest stop mass mi^ = (100,125,130,135,150,160) GeV 
and showing the limits in the Ot-mi^ plane. The individual constraints from the 77 , Z 7 , ZZ, 
and WW channel are represented (no bounds from hh channel can be obtained from the 
current data). For comparison we also display the region excluded by indirect constraints 
coming from Higgs signal strength measurements, electroweak precision data, and vacuum 
stability (see Appendix]^. In contrast to the previous section in which we defined ideal or 
benchmark branching ratios, here we account for the realistic branching ratios, computed 


at each point in the physical stop parameter space as described in Sec. IIG 
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FIG. 8 : Excluded regions in the physical stop parameter space. Here we show the limits from 
diboson resonance searches in the plane for fixed = (100,125,130,135,150,160) GeV, 

including 77 (blue hatched), Zj —)> £+£“7 (green), ZZ —^ M (red), WW —)• 2i2i' (orange). We 
also display the union of indirect constraints (gray) shown in Fig. 11 The strongest limit from 
ATLAS or GMS in each channel are shown. Finally, we indicate the stop parameters for which 
= 0 (red dashed contour). 
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Assuming the stoponium forms and annihilate decays, we observe that light stops with 
masses mi^ < 130 GeV are excluded by the combination of diboson resonance searches 
and the indirect constraints. Stoponium searches alone probe stops lighter than about 125 
GeV. Remarkably, a strong complementarity is seen between the 77 resonance searches and 
the indirect constraints (particularly the Higgs coupling constraints), with the 77 channel 
probing the remaining open stop blind spot region corresponding to small or vanishing 
coupling Moreover, we also see that primarily left-handed stops with masses < 

170 GeV are constrained by Z 7 resonance searches. 

Again, we emphasize that the limits we present here rely on the A^ = 300 MeV pa¬ 
rameterization of the stoponium wavefunction at the origin from Ref. [23] and the inclusion 
of the hrst two S'-wave states PDl. While it is likely that this estimate is conservative, 
particularly in light of the possible additional signal coming from the higher excited states, 
it is worthwhile to examine how the limits change if different assumptions are made. For 
example, if we instead include only the contribution of the IS state, the expected limit on 
the stoponium mass will be weaker by about 20 GeV, as can be seen by examining Figs. 
and Alternatively, if one assumes the pure Goulomb potential, the expected limit on the 
stoponium mass will be stronger by about 70 GeV. In any case, it is clear that stoponium is 
now probing the hypothesis of light stops beyond the LEP limits and will explore uncharted 
territory with the next run of the LHG. 


V. EXAMPLE SCENARIOS OF INTEREST FOR STOPONIUM 

We now discuss two motivated examples scenarios in which stoponium provides a comple¬ 
mentary probe to direct stop searches and indirect tests. The condition ([^ will generically 
hold when the stop does not possess an unsupressed 2-body decay. This can happen due 
to kinematics or small couplings, and to demonstrate this we will focus two examples - 1 ) 
a “canonical” case of R-parity conservation with a bino LSP and a stop NLSP, and 2) a 
R-partiy violating case with a stop LSP which decays to a pair of jets through the UDD 
operator. 


A. R-Parity conserving bino LSP and stop NLSP 

A canonical benchmark scenario with light stops consists of a stop NLSP and a neutralino 
LSP, under the assumption that R-parity is conserved. In this case the neutralino is stable, 
leading to the signature of missing energy at colliders. As is well known, in this scenario 
the stop width is very sensitive to the mass spectrum, and several kinematic regions can be 
distinguished: 
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1. > mt + m^o: The stop decays via the 2-body process ii —)■ 

II. mw + mb + < mi^ < mt + : The stop decays via the 3-body process ti —)■ Whx^■ 

III. rrib + m^ < mi^ < miyThe stop decays via the 4-body process ii —)■ bff'x^, 
or possibly through the flavor-violating 2-body decay ii —>■ cy^. 

To illustrate these features, in the left panel of Fig. |^we have plotted the decay width of 
the stop as a function of its mass for the case of a 10 GeV bino LSP (see the Appendix 
for the 3-body partial decay width of the stop). For comparison, in the same plot we show 
the stoponium binding energy and annihilation decay width. As the stop mass approaches 
the boundary between 2- and 3-body decays, we observe that F^ becomes larger than the 
natural stop width. Therefore, in region I above, the stop generally decays too quickly 
for stoponium resonance signatures to be visible, although it may be possible near the 
edge of region I where phase space suppression is relevant. However, in regions II and III 
the dominant stop decay channels are 3- and 4-body (or suppressed 2-body), respectively, 
allowing both sufficient time for the formation of the stoponium and the dominance of the 
bound state annihilation decays. 

Therefore, regions II and III are prime targets for resonance searches for stoponium, which 
is clearly illustrated by the “phase” diagram in the stop - neutralino mass plane shown in the 
right panel of Fig. We represent the current exclusion limit from stoponium resonances by 
the blue hatched region; this region is excluded for the given lightest stop mass for all other 
stop parameters (see Fig. |^. For comparison, we have overlaid the current ATLAS limits 
from direct stop searches [6], represented by the orange dotted line. While these direct limits 
naively appear to cover most of light stop region, it is important to emphasize that they rely 
on the assumption of a 100% branching ratio to the hnal state under consideration. If this 
assumption is relaxed, the limits are weakened considerably. For example, in region III if 
the stop has comparable branching ratios in the 4-body bff'x^ and 2-body cx^ channels, a 
signihcant portion of parameter space opens up [2H] , and only the orange-shaded region can 
be excluded regardless of the assumption on the branching ratio. In contrast, constraints 
from stoponium resonances are already starting to unambiguously probe this region. 

It is also worthwhile to note that the “stealth” stop region - the phase space-suppressed 
part of region I near ~ rrit, m^o > 0 - is also potentially amenable to stoponium searches 
(within the red-shaded region in the bottom part of Fig. This region is challenging to probe 
with direct stop searches, and several dedicated observables such as spin correlations [2^130] . 
top pair production rates [HU 132] and special kinematic variables |33| have been studied. 
However, they primarily rely on precision measurements and calculations, which are often 
subject to subtle uncertainties. On the other hand, the stoponium resonance searches are 
cleaner and less ambiguous in most cases. Thus, the stoponium can provide an important 
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FIG. 9: R-parity conserving bino LSP and stop NLSP: (Left): Comparison of the stop decay 
width (red), stoponium annihilation rate (blue) and stoponium binding energy (green). Here we 
have assumed a pure bino LSP with Mi = 10 GeV, = 1 TeV, and 9t = 0.18. The vertical line 
separates the two- and three-body decays of the stop. (Right): In the red shaded region above 
each red line, the stoponium annihilation is quicker than the individual stop decay, and searches 
for stoponium resonances are possible. Three benchmark scenarios are presented as red lines: 
the top two have almost ideal 77 branching ratios (defined in Sec. Ill A) with 9t ~ 1.36 (top), 
0.18 (middle), and the bottom one with 9t = 0.12 (bottom). For comparison, we also display a 
union of ATLAS limits on the stops (orange dotted), which assume 100% branching ratios to the 
final states under consideration [ 6 j as well as a branching ratio-independent exclusion limit taken 
from Ref. [28] (orange shaded). The blue-hatched region is excluded from stoponium Fig. ^ 


complementary probe of the stealth region. In this region, stops with a mostly left-handed 
mixture (small 9t) can be probed via stoponium over a wider range of stop masses (see Fig. 
right) because the individual stop decay width is smaller with a smaller 6 t. In any case, 
stoponium does not yet have sensitivity to stealth stops (the blue-hatched region extends to 
about 130 GeV), but can potentially be sensitive to a significant portion of this challenging 
region in the early part of Run 2 with C)(10)/fb of data (see Fig. [^. 

Let us finally comment on the case in which the LSP is a very weakly interacting particle 
such as the gravitino. In this scenario, the stop NLSP is typically long lived due to its weak 
coupling to the gravitino, and therefore stoponium annihilation decays can be relevant. 
For instance, taking a low supersymmetry breaking scale y/F = 10 TeV (corresponding 
to an essentially massless gravitino), a 200 GeV stop has a decay width T{ti —)■ tG) ~ 
10“® GeV [3l|, which is much smaller than the typical stoponium annihilation width. In 
practice, however, this scenario is already strongly constrained by direct searches covering 
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OTjj [GeV] 

FIG. 10; R-parity violating stop LSP: We display contours for two values of the stop width, 
= 1 GeV (red) and F^^ = 10“^ GeV (blue) for typical values of stoponium binding energies and 
annihilation decay widths. We also display a contour of the stop lifetime = 0.1 mm. Stops 
with lifetimes in this range and longer are subject to stringent LHC searches for long-lived charged 
particles [351111]. Stops with prompt RPV decays in the red-shaded region are also constrained 
from RPV searches |36|. Stoponium can therefore be relevant for a broad range of RPV couplings, 
a” < 10 “^, and the current exclusion from stoponium resonances is shown as blue-hatched. 

both prompt and displaced decays [35|, and only a small window aronnd the stealth stop 
parameter region for prompt decays is still open. The stoponinm searches will be able to £11 
this open region. 


B. RPV stop LSP 

Another benchmark scenario of interest is a stop LSP that decays to a pair of jets dne 
to a small UDD R-parity violating conpling, which we cnstomarily denote as A". Snch a 
RPV stop is challenging to probe with direct searches at LHC. Despite its large prodnction 
cross section, the signatnre of fonr jets (paired dijet resonances) is difficnlt to disentangle 
from the QCD backgronnd, and stops as light as ~ 100 GeV are still allowed HHiiaiMI 
althongh some range of the stop masses 200 — 350 GeV has recently been exclnded |37|. 

Stoponinm resonance searches provide another means of probing this interesting scenario. 
Since the decay of the stop is throngh a 2-body process, the RPV conpling mnst be somewhat 
snppressed in order for the condition ([^ to apply. In fact, this may be motivated in explicit 
models of R-parity violation which aim to protect the lifetime of the proton throngh other 
symmetries. For example, in the scenario of MFV SUSY jSHlEn], the dominant RPV conpling 
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is \sb ~ 


0(10“^), which is small enough for ^ to hold. Note also that if the RPV coupling 
is extremely small, A" < 0(10“^) then the stop will have a macroscopic lifetime. In this 
case, there are additional strong constraints from searches for displaced dijets [35l 00] or 


long-lived colored/charged particles [3510T]. In Fig. 10 we display the RPV stop LSP phase 
diagram. We conclude that stoponium resonance searches can provide new complementary 
probe of this scenario for a broad range of couplings, 10“^ < A” < 10“^. 


VI. OUTLOOK 

A natural supersymmetric solution to the hierarchy problem predicts light stops, which 
can lead to a variety of novel phenomena. Besides the suite of direct searches at the LHC 
(which depend on the SUSY spectrum) and indirect tests such as Higgs couplings, preci¬ 
sion electroweak data and vacuum stability requirements, it is possible that light stops can 
manifest through a stoponium annihilation decay leading to a resonance signature in the 
77 , Zy, ZZ, lUlU and hh channels. This will happen if the stop does not have an unsup¬ 
pressed two body decay and naturally occurs in several motivated scenarios. 

In this work we have derived new limits on light stops from ATLAS and CMS diboson 
resonance searches using Run 1 data. Our limits are derived using the empirical quarkonia 
potential model of [21] and assuming the hrst two S'-wave states contributing to the res¬ 
onance signal, which we have argued is conservative. Assuming the stoponium can form 
and annihilate decay, light stops below about 130 GeV are excluded by the combination of 
constraints on stoponium resonances and other indirect constraints, such as Higgs couplings 
measurements. Notably, we have demonstrated that the 77 channel is especially complemen¬ 
tary to these indirect probes, as the resonance signature in this channel is maximized in the 
blind spot region where the Higgs-stop-stop trilinear coupling is small and these constraints 
are weakest. Furthermore, Z 7 resonance searches provide even stronger limits for primarily 
left handed stops, limiting stop masses below about 170 GeV. In the long term, the LHG 
experiments will be able to probe stops in the 300-400 GeV range with a high-luminosity 
run with 3/ab at 14 TeV. 

We have also discussed the implications of these searches for two specihc SUSY scenarios 
with lights stops. In the case of a bino LSP and stop NLSP, stoponium annihilation decays 
can be relevant since the stop width is naturally suppressed over a wide range of parameters. 
Stoponium resonance searches provide a relatively clean and unambiguous probe, and can 
give a handle on the challenging compressed and stealth stop regions of parameter space. 
Instead for the case of a R-parity violating stop LSP decaying to two jets, stoponium already 
provides unique sensitivity to stops lighter than about 130 GeV and RPV couplings in 
the range 10“^ < A" < 10“^, which is difficult to probe directly due to the large QGD 
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backgrounds. 

Looking forward, it is clearly of interest on the theoretical side to bring under better 
control the uncertainties coming from our imprecise understanding of the stoponium bound 
state system. A modern study of empirical quarkonia potential models along the lines of 
Ref. [21] is warranted, and perhaps lattice studies could provide further insight into the 
non-perturbative matrix elements entering into stoponium production and decays. A more 
detailed investigation into the decay patterns of the excited states would also be helpful 
in order to understand their contribution to the resonance signals and perhaps uncover 
additional signatures of stoponium. On the experimental side, stoponium clearly provides 
a well motivated target for resonance searches, which should be a high priority during the 
next run of the LHC. 
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Appendix A: Stop and sbottom sector conventions 

In terms of the gauge eigenstates {tL, tR), the stop mass matrix is given at tree level by 

I rriQ^ +m^ + rritXt 

y mtXt 

where are the left- and right-handed squark soft mass parameters, Xt = At — 

/i/tan/?, with At the soft trilinear coupling, /r the supersymmetric Higgs mass parameter, 
and tan/? the ratio of up and down type Higgs vacuum expectation values, and D\ = 
m| cos2/? (I — |s^), = m| cos2/? Eor simplicity, we will assume all parameters 

are real. The physical stop mass eigenstates are related to the gauge eigenstates through 
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the orthogonal transformation: 



cos 9t — sin 9t 
sin 9t cos 9t 



(A. 2 ) 


where the mixing angle 9t satishes 


tan 20 + = 


2mtXt 


t ’ 




(A.3) 


which we take in the range [ 0 , vr/ 2 ]. 

The sbottom sector is described in a similar fashion. The sbottom mass matrix can 


be obtained from Eq. lA.ll with the replacements —>■ rrib, mu^ —)■ Xt ^ = 


Ab — /itan/3, cos 2/3 (—| + and —)■ = m| cos 2/3 (— 


Appendix B: Indirect constraints on light stops 


In this appendix we summarize the existing indirect constraints on light stops coming 
from Higgs signal strength measurements, precision electroweak data, and vacuum stability. 
In the MSSM, there will be additional constraints on the stop parameters coming from the 
requirement of obtaining a 125 GeV Higgs mass, but we will remain open to new contribu¬ 
tions to the Higgs quartic coming from physics beyond the MSSM. Furthermore, if Higgsinos 
or gluinos are light there can be further constraints from flavor physics, which we will not 


include here. Several summary plots of these constraints are shown in Fig. 11 


1. Higgs signal strength data 

Light stops modify the hgg and hyy couplings at one loop and are therefore subject to 
constraints by ATLAS and CMS measurements of the Higgs signal strength parameters BZI- 
To derive the exclusions, we use the method of Ref. [13], with the recent updates in Ref. P], 
which include the latest Run 1 results as of summer 2014. In this approach, the raw signal 
strength data are combined to derive a function for the 77 , VV, and h/r channels as 
a function of two combined production mode signal strengths: 1 ) gluon-gluon fusion and 
tth (ggF-|-ttH), and 2) vector boson fusion and associated production (VBF-I-VH). This 
yields eight composite signal strengths, In terms of the ratios r, = r(h —)■ 
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i)/V{h —>■ i)sM, the signal strength predictions are given by 


■'ggH+ttH ^ ^gg ^77 -VBF+VH 

“77 

,7ggH+ttH _ -ggH+ttH _ -ggH+ttH _ 
h^VV r^hh A^rr 


^77 

rr 


-VBF+VH _ -VBF+VH _ -VBF+VH _ 
H'VV f^bb h'rr 


(B.l) 


where we have dehned rr = ~ 1 + (vgg — l)BR(h —)■ gg)sM with BR(h —)■ gg)sM — 

0.085. 

One loop sbottom and stop exchange lead to the following expressions for Vgg, r^y^. 


^gg = 


1 + 


^/2vC{r) 

^SM 




^hOhq 2 ' + 
"h'l 
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where C'(r) = 1/2, (i(r) = 3, Qi = 2/3, = —1/3, and 1.3, 6.6. The 

Higgs-squark trilinear conplings Aj can be fonnd in, e.g., Ref. |1S]. The scalar loop fnnction 
24 o(^) is dehned in Ref. [10] and approaches 24o(m) —>■ 1/3 in the limit m 3> rribp. 

With these ingredients we constrnct the fnnction. As Vgg and are correlated, we 
dehne a one parameter best-ht region in the physical stop parameter space at 2a C.L. by 
demanding < 4. 


2. Electroweak precision data 

Regarding the precision electroweak data, the largest effect of the stops and sbottoms is 
to provide a new one-loop contribntion to the p parameter m- 


Ap = 


3Gf 


8 V 2 


TT^ 


StCtFoim 




) + c/Fo(m?,m?J + Fo{m 




(B.4) 
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where Gp^ = 2\/2v'^ and we have defined the function 


Fo(x, y} = x + y -log 

x-y y 


(B.5) 


We have used the left-handed sbottom mass, in Eq. (B.8) with the zero sbottom mixing. 
We apply the constraint from Ref. 


Ap= (4.2 ±2.7) X 10"h 


(B.6) 


3. Vacuum stability 


In the field space of { ~ h, f/j}, we study the (meta)stability conditions for 

our electroweak vacuum characterized by {h)/ a/2 = v = 174 GeV and vanishing stop 
vevs [4^ [50] . The tree-level potential that we use is (in terms of real scalar fields with 
1/v^ normalizations factored out) 


T/O ^ 2 1,2 I ^ 2 22 

V = -mirh + -m ^ 
2 ^ 2 
J2 


^ ([iLtRf ± {htLf ± {htR)‘ 


Vt ^ 


‘'Q 


32 V 3 ^ 




+ 

32 V 


72 


+ 24 


^R 


ytXt ~ ~ , 4A, 4 

± — - 1 =^- htijtR ± h.c. ± . 

a/2 4 


(B.7) 


We have assumed that the effects of are negligible, appropriate in the large regime, and 
have substituted sy —)■ 1, C /3 —0, and —)■ h/\/2. To have a proper electroweak vacuum 

and the measured Higgs mass, we set rrijj = rrih = 125 GeV, and 6X = ^ 

In particular, we assume the presence of some physics, either radiative stop corrections in 
the MSSM or physics beyond the MSSM, that generates an appropriate correction to the 
Higgs quartic coupling, 6X. To arrive at the above form of the gs terms, we assume that the 
stops are aligned or anti-aligned in the SU{3)c color space. 

The negative Xt - linear term can induce new minima when all three scalar fields obtain 
vevs. We numerically scan the field space to determine all local minima, and we use the 
CosmoTransitions program [51] to calculate the tunnelling rate between our electroweak 
minimum and any deeper charge-color breaking minima. The classical action for this tunnel¬ 
ing, S'e > 400 is required for the metastability of our vacuum. We find that the parameter 
space with large values of the Higgs-stop-stop coupling, is constrained from the 

requirement of vacuum (meta)stability. 
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Our vacuum stability bounds rely on the simplifying assumptions discussed above, and it 
is possible that stronger constraints can be obtained if some of these assumptions are relaxed 
(e.g. Hd vevs or misaligned stop vevs), or when one accounts for one-loop corrections to the 
potential and thermal effects (for a useful comparative study, see Ref. |50]). Nevertheless 
we believe our bounds give a qualitatively correct representation of the current limits, and 
we observe that for light stops the vacuum stability limits are generally weaker than those 
coming from Higgs couplings. 

When the left-handed stop is light, we also check whether the mass of the accompanying 
left-handed sbottom is positive (and evading LEP bounds). For the given stop parameters 
and heavy enough right-handed sbottoms, the left-handed sbottom is the heaviest when the 
sbottom mixing is zero (the equality below holds in this case) 

4- - ml - C2p m^. (B. 8 ) 

Throughout in this work, we use this left-handed sbottom mass with the equality sign by 
assuming the zero sbottom mixing. 


Appendix C: Stop 3-body decay: mt + m^a > mr > mw + mi, + m^o 


The partial decay width for t ^ t*x^ ^ Wbx^ can be written as 


ii^Wbx° 


J {mw 


ds 


(5 _ miy + mlTl (C.l) 


where we have dehned the “off-shell” partial widths: 
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(C.2) 


Here, we have dehned the couplings. 


Sl = NliS, + y,Nl,c,, 

1 '/v W 
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9 ^12 Ct + yt Ni 4 St, 
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FIG. 11: Indirect constraints on light stops: Here we display for light stop masses = 
(100,130,160,200,300) GeV the constraints in the 9t - mi^ plane coming from Higgs signal strength 
measurements (gray), precision electroweak data (orange), vacuum stability (yellow), tachyonic 
sbottoms (blue) and sbottoms below the LEP kinematic reach (light blue). The white regions are 
currently unconstrained. The red dashed line indicates where the Higgs-stop-stop coupling 
vanishes. Here we have hxed tan/3 = 10 and the sbottom mixing to zero {X^ = 0). 
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with Nij the neutralino mixing matrix elements, and the kinematic function X{a,b,c) = 
— 2ab — 2ac — 2bc. Notice that the three-body partial decay width reduces to the 
two-body one in the appropriate kinematical regime through the use of the narrow-width 
approximation. 


Appendix D: Limit Extrapolation to Future Searches 


To estimate the future prospects for probing stoponium at the LHC we extrapolate the 
current expected limits of the diboson resonance searches to ^/s = 14 TeV with integrated 
luminosities C of 30/fb and 3/ab. To accomplish this, we employ several simplifying as¬ 
sumptions. For the uncertainties we consider two extreme cases: (1) statistical uncertainties 
dominate, and (2) systematic uncertainties dominate and improve in proportion to a/Z. The 


results in the main text of Section are obtained with the former assumption, but we will 
offer a comparison between the two cases below. 

Statistical uncertainty dominant: At a specihed collision energy and integrated luminos¬ 
ity denoted collectively by a superscript i, we have 


(T* 1 1 

oc , (D.l) 

where cr^ound upper bound on the signal cross-section, cr^ is the background cross- 

section, and is the number of background events after all cuts. The latter two quantities 
are related as A’|j = = aBVqqC‘e\ where cfb is the partonic background cross-section 

(which does not depend on i), O is the integrated luminosity, e* is the selection efficiency, 
and is the parton luminosity for initial partons ab, dehned as 

ru= dx-Ux,Q^)MT'/x,Q^), r‘ = ya. (D.2) 

Jt' ^ ^ 

we have taken Vg^ since the dominant backgrounds in the most sensitive chan¬ 
nels, 77 and ZZ, arise from qq initial states. In our numerics we take the up quark parton 
luminosity for simplicity, as it will dominate in pp collisions. 

We furthermore assume a constant efficiency, e* = e, which is reasonable since the signal- 
to-background ratio near a resonance is largely determined by the resonance mass. This is 
similar to gluino pair searches |52] and perhaps more generally useful [53] . 

Taking ratios, we can extrapolate the current bounds to those for other collision energies 


In Eq. (D.l 
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and luminosities in a simple way in terms of parton luminosity ratio and luminosity ratio 
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Again taking ratios, we are able to extrapolate the current limits: 
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Comparing Eq. (D.3) and Eq. (D.5), we see that if errors are dominantly from systematics, 


one obtains a weaker cross-section bound by a factor yJVlq/V/q. Numerically, this equals to 
1.45 and 1.73 for = 160, 400 GeV at 14 TeV collision for q = u. These results roughly 
show how sensitively the future bounds may depend on the assumptions on errors. 


Appendix E: Resonance Searches 

In Table [I], we collect the latest experimental results on diboson resonance searches. The 
searches shown in bold give the strongest constraints in each diboson channel and are used 
to derive limits on stoponium. The CMS and ATLAS 77 limits are similar but constrain 
complementary mass ranges due to different statistical fluctuations, and we therefore show 
both results. The CMS Z 7 is found to be somewhat stronger than that of ATLAS, but due 
to the importance of this channel, we show both results and encourage more careful analysis. 
All searches use 8 TeV datasets except for the ones with ~24/fb and ~10/fb, in which case 
the combined 7-1-8 TeV dataset is used. The CMS ZZ —)■ 4£ result is stronger than the 
ATLAS one partly because CMS presents combined results of all production channels and 
uses the 7-1-8 TeV dataset. The WW —)■ ^v2j result has similar sensitivity to that of the 
WW —)■ 2^2v, but we show only the latter result for simplicity. The searches listed in the 
last four lines give weaker constraints. Lastly, ATLAS 77 and Z 7 present limits on the 
fiducial cross-section, so we assume and unfold a constant hducial efficiency for each channel 
as mentioned in the relevant part of text. 
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channel 

CMS 

ATLAS 

77 

19.7/fb 

HIG-14-006 [M] 

20.3/fb 

1407.6583 [55] 



19.7/fb 

HIG-14-031 [SB] 

20.3/fb 

1407.8150 [ST] 


zz 

24.8/fb 

1312.5353 [SB] 

20.7/fb 

CONF-2013-013 

PI 

WW iiy2j 

24.3/fb 

HIG-13-027 [60] 


- 


WW 2i2u 

24.3/fb 

1312.1129 [SI] 

20.7/fb 

CONF-2013-067 

ra 

h/i —J- 6677 

19.7/fb 

HIG-13-032 [63] 

20.0/fb 

1406.5053 [SI 


ZZ 2i2j 

10.4/fb 

1304.0213 [65] 

20 .0/fb 

CONF-2014-039 

[SS] 

ZZ 2i2u 

10.4/fb 

1304.0213 [SS] 


- 


hh —)■ Ah 

17.9/fb 

1503.04114 [ST] 

19.5/fb 

CONF-2014-005 

PI 

hh WWbb 


- 

20.3/fb 

1312.1956 [S9] 



TABLE I: The latest data for heavy resonance searches used in this paper. Those in bold are used 
to derive limits on stoponium. Other results give weaker constraints. 
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